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Abstract 

With the HIJING/BB v2.0 heavy ion event generator, we explore the phenomenological con- 
sequences of several high parton density dynamical effects predicted in central Pb+Pb collisions 
at the Large Hadron Collider (LHC) energies. These include (1) jet quenching due to parton en- 
ergy loss {dE/dx), (2) strangeness and hyperon enhancement due to strong longitudinal color field 
(SCF), and (3) enhancement of baryon-to-meson ratios due to baryon-anti-baryon junctions (J J) 
loops and SCF effects. The saturation/minijet cutoff scale pq{s,A) and effective string tension 
k{s, A) are constrained by our previous analysis of LHC p + p data and recent data on the charged 
multiplicity for Pb+Pb collisions reported by the ALICE collaboration. We predict the hadron 
flavor dependence (mesons and baryons) of the nuclear modification factor i?AA(PT) and emphasize 
the possibility that the baryon anomaly could persist at the LHC up to px ~ 10 GeV, well beyond 
the range observed in central Au+Au collisions at RHIC energies. 

PACS numbers: 12.38.Mh, 24.85. +p, 25.75.-q 
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I. INTRODUCTION 



With the commissioning of the Large Hadron CoUider (LHC), it is now possible to test 
dynamical models of multiparticle production in nuclear collisions up to a center-of-mass en- 
ergy (cm.) per nucleon ^snn = 5.5 TeV. Charged particle densities at mid-pseudorapidity, 
{2dNch/ drj) / Npi^j-t, and their dependence on energy and centrality are important for under- 
standing the mechanism of hadron production and especially the interplay of soft (string 
fragmentation) and hard (perturbative quantum chromodynamics) scattering contributions 
at an order of magnitude higher energies than were extensively studied up to now at the Rela- 
tivistic Heavy Ion Collider (RHIC/BNL). The rate of parton-parton and multi parton-parton 
scattering and soft processes are strongly correlated with the observed particle multiplicity 
(related to the initial entropy and the initial energy density generated in nuclear collisions 
over spacetime volumes up to 10^ greater than in elementary p + p collisions). 

Using the constraints on dynamical parameters from data at RHIC energies, there remain 
relatively large uncertainties (~ a factor of two) on predictions for charged hadron multi- 
plicities at the LHC energies, see Refs. [l-16|. The first data on inclusive charged particle 



distributions from the LHC in Pb+Pb collisions are now available [17H22|. Recently, AT- 
LAS analyzed event samples from central (0-10%) Pb+Pb collisions at ^/s^ = 2.76 TeV 
in which the transverse energy of the most energetic jet (leading jet) E^^ > 100 GeV and 
where a recoiling jet with transverse energy E^^ > 25 GeV could be found at azimuthal 
separation A$ > 7r/2. In comparison with pp collisions, ATLAS observed (within a narrow 
cone) in Pb+Pb collisions a strong reduction of the number of recoiling jets carrying a large 
fraction of the maximum available jet energy x = Et^/ Et^ > 0.6 21[. This is consistent with 
the expected strong quenching of jets due to partonic energy loss in ultradense quark gluon 
plasma (QGP). In addition to the tomographic interest in jet quenching high px observables, 
there is considerable interest on the impact of mini-jet quenching as an additional final state 
interaction source of the observed bulk multiplicity /entropy production. The LHC data 
on nucleus-nucleus {A + A) collisions may lead to an improved theoretical understanding 
of ultra-dense multiparton plasma based on a quantum chromodynamics (QCD) approach 



0,0, 



23 



25|. 



To isolate quark-gluon plasma signatures in heavy-ion collision data, it is imperative to 
work in a framework that can simultaneously account for nucleon-nucleon collisions {p + p) 
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in the same energy range. In fact p + p reactions can not be considered as elementary at the 



LHC and their study has aheady revealed possible interesting new physics |26|. We utilize 
the HIJING/BB2.0 heavy-ion generator as a phenomenological tool that can systematically 
relate observables in p+p, p+A, and A + A collisions in terms of a common framework based 
on a two-component picture of the dynamics involving soft longitudinal beam jets, multiple 
minijets, and rare hard pQCD processes. The Monte Carlo code inherits the phenomeno- 



logical successes of the LUND 



27( 1 and dual parton model (DPM) |28| string excitation and 



hadronization mechanisms as well as the hard processes encoded in the PYTHIA 29| event 
generator. Multiple low p^ < Qs{x,A) transverse momentum color exchanges excite the 
incident baryons into longitudinal strings that fragment due to color neutralization into an 
array of physical hadron resonance states. Hard processes are modeled as high px kinks 
in the strings and hadronize via the well-tested LUND scheme in e+e~ and ep processes. 
The HIJING model 30| uses a variant of the string excitation LUND and DPM models 



constrained to lower energy pp data. 

In heavy ion collisions, the novel "nuclear physics" is due to the nuclear modification of the 
parton distribution functions, the possibility of multiple longitudinal flux tube overlapping 
leading to strong longitudinal color field (SCF) effects refered to as color ropes or glasma. 
Strong fields also lead to enhanced strangeness production [3l|. The effect is especially 
enhanced for production of baryons with multiple strange quarks as discovered first at the 
SPS/CERN |32|. Recently, the data on the production of strange particles with one, two or 
more strange quarks for p + p collisions at the Tevatron j33| and for p + p collisions at the 



LHC 
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35| have been reported. In Ref. 33[ it was shown that the ratio S /A and Q /A 



rise at low pr and the ratio reaches an unexpected plateau at > 4 GeV/c, which persist 
up to rather large 10 GeV/c. This result is a major motivation for the present work. 

We address these and other issues within the framework of the HIJING/BB v2.0 model 



36-l38|. A systematic comparison with data on p+p and p+p collisions in a wide energy range 
36l | revealed that minijet production and fragmentation as implemented in the HIJING/BB 
v2.0 model provide a simultaneous and consistent explanation of several effects: the inclusive 
spectra at moderate transverse momentum [px < 5 GeV/c), the energy dependence of 
the central rapidity density, the (strange)baryon-meson ratios and the baryon-antibaryon 
asymmetry. Specifically, in this paper we extend our studies to the dynamic consequences of 
initial and final state saturation phenomena, gluonic JJ loops, and jet quenching to particle 
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production in central (0-5%) Pb+Pb collisions at the LHC Our study aims to investigate 
a set of observables sensitive to the dynamics of the collisions, covering both longitudinal 
and transverse degree of freedom. We present here the pseudorapidity multiplicity density 
per participant pair, and their centrality and energy dependence. Our study reveals a 
possible violafou at ^ 2.76 TeV of the „„otee.ved in „.uMpa.tic.a 

production of charged particles up to top RHIC energy [39|, |40|. Predictions for the nuclear 
modification factors (NMF) of charged particle (-Raa) and of identified particles (-Raa) 
also discussed. We emphasize the role of SCF effects and gluonic JJ loops in understanding 
the meson/haryon anomaly at the LHC energy, which manifest in the (strange)baryon-meson 
ratios at intermediate and large transverse momenta and in NMF of identified particles 

(^aa)- 



II. OUTLINE OF HIJING/BB V2.0 MODEL. 



A detailed description of the model can be found in Refs. 36l438l|. The HIJING model 
is based on a two component geometrical model of minijet production and soft interaction 
and has incorporated nuclear effects such as shadowing (nuclear modification of the parton 
distribution functions) and quenching, via final state jet medium interaction. In the 



HIJING/BB v2.0 model 36|, |37| we introduce new dynamical effects associated with long- 



range coherent fields (i.e., strong longitudinal color fields, SCFs), including baryon junctions 



and loops [38|, |4]| . 



Saturation physics is based on the observation that small-x hadronic and nuclear wave 



functions, and thus the scattering cross sections as wel 
momentum scale known as the saturation scale (Qsat) 



, are described by the same internal 



to LHC 7 TeV 



42|. A recent analysis of pp data up 



43| has shown that, with the factorized gluon fusion approximation the 



growth of the dNch/drj can be accounted for if the saturation scale grows with center-of-mass 
system (c.m.s.) energy as 



sat V 



Qlis/so) 



\CGC 



(1) 



with AcGC ~ 0.115. The saturation scale is also increasing with atomic number as A^^^ 44 1 
and it was argued that the effective string tension (k) of color ropes should scale with Qs^t 



44 



45|. 



However, in HIJING the string/ rope fragmentation is the only soft source of multiparticle 
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production and multiple minijets provide a semi-hard additional source that is computable 
within coUinear factorized standard pQCD with initial and final radiation (DGLAP evolu- 



tion [46|). In order to achieve a quantitative description, within our HIJING/BB framework 
we will show that combined effects of hard and soft sources of multiparticle production can 
reproduce the available data in the range 0.02 < ^/s < 20 TeV only with a reduced depen- 
dence of the effective string tension on ^/s. We find that the data can be well reproduced 
taking 

k{s) = Ko (s/so)'-"' GeV/fm ^ Qo. g,at(s), (2) 

where kq = 1 GeV/fm is the vacuum string tension value, Sq = 1 GeV^ is a scale factor, 
and Qok = i^o/Qo is adjusted to give k = 1.88 GeV/fm at the RHIC energy ^/s = 0.2 TeV. 



Our phenomenological k,{s) is compared to Qsatl-^) Fig. 1 from Ref . |36l| , where k = 1.40 
GeV/fm at = 0.017 TeV increases to k = 3.14 GeV/fm at = 14 TeV. 

.„ a strong longitudinal color electric field, the heavier flaw suppression factor varies 
with string tension via the well-known Schwinger formula |47J, 

r„o expf-"<^^«-'"''l<l (3) 



^qq \ f^O 

for Q = qq, s, c or 6 and q = u, d. Therefore, the above formula implies a suppression of 
heavier quark production according ton:(i:qq:s:c~l:l: 0.02 : 0.3 : 10~^^ for the 
vacuum string tension Kq = 1 GeV/fm. For a color rope , on the other hand, the average 
string tension value increases due to in-medium effect. With these increase string tensions 
the QQ flavor pair production suppression factors, 7qq, approach unity in A + A collisions. 
In our model this is the main mechanism for strange meson and hyperon enhancement. 
This increase is quantifled in our calculations at RHIC and LHC energies, using a power-law 
expression, 

<s,A)j,Hic = ^{s)A'-''' (4) 
n{s,A)^HC = ^{s)A'-''' (5) 

where the exponent has been flxed to provide a good description of the charged particle 
densities at mid-pseudorapidity (see below). 

3 

Moreover, the measured charged multiphcity grows faster (~ y/s ' ) in nucleus-nucleus 



collisions than it does for protons (~ 
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The energy dependence of the multiplicity 
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is an experimental probe of possible x dependence of the saturation scale Ql^^ |49|. Also, 
the energy density for most central (0-5%) collisions is larger by about a factor of 3 at LHC 



than the corresponding one at RHIC 48|. Since the energy densities are computed from the 



square of the field components 



50| they are proportional with n[s, A^. These experimental 



facts indicate that we have to use different exponent at RHIC and at LHC. The problem 

1 

was also addressed in CGC models, where ~ Qsp(^^ )^''^ • smaller x (corresponding 
to LHC energy) l/7cs is larger, and thus the nuclear enhancement of Q^^^ is expected to be 



larger than at top RHIC energy [49|]. 



It was shown recently 



5l| that charged-particle multiplicities measured in central heavy 



ion collisions at high energies may not directly determine the initial conditions as predicted 
by CGC models. In its simplest implementation it was generally assumed that there is 
a direct correspondence between the number of partons in the initial state and the num- 



ber of particles in the final state 
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431 ]. However, final-state nonequilibrium production 



mechanisms generally increase the final multiplicity. In fact, in recent HIJING2.0 (which 
includes modern structure functions but without SCF and JJ loops) the enhancement of the 
multiplicity due to quenching of final minijets had to be turned off not to over predict the 
ALICE dNch/dy j^. In HIJING/BB2.0 used here with supersaturated Duke-Owens parton 
distribution functions we find a consistent description of both p + p and Pb+Pb multiplicity 
when taking into account the extra multiplicity generated by minijet quenching. 

As mentioned, in HIJING the string/rope fragmentation is the only soft source of multi- 
particle production and multiple minijets provide a semihard additional source that is com- 
putable within coUinear factorized standard pQCD with initial and final radiation (DGLAP 



evolution 



46|). Within the HIJING model, one assumes that nucleon-nucleon collisions at 



high energy can be divided into soft and hard processes with at least one pair of jet pro- 
duction with transverse momentum, px > Po- A cutoff (or saturation) scale po{s,A) in the 
final jet production has to be introduced below which the high density of initial interactions 
leads to a nonperturbative mechanisms which in the HIJING framework is characterized by 
a finite soft parton cross section cXsoft- The inclusive jet cross section cijet at leading order 
(LO) I52I is 



fs/4 -j^ 

ajet = I dp^dyidy2 
hi 



dcr jet 



2 dpl,dyidy2 ' 



(6) 
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where, 



- K 2_^Xifa{Xl,PT)X2fb{X2,PT)- 



dpldyidy2 



a.b 



dt 



(7) 



depends on the parton-parton cross section a"''' and parton distribution functions (PDF), 
fa{x,p\). The summation runs over all parton species; yi and y2 are the rapidities of the 
scattered partons; Xi and X2 are the light-cone momentum fractions carried by the initial 
partons. The multiplicative K factor [K 1.5 — 2) accounts for the next-to- leading order 
(NLO) corrections to the leading order (LO) jet cross section 53|,|5J|. The parton distribu- 



tions per nucleon in a nucleus (with atomic number A and charge number Z), fa/A{x, Q^, r), 
are assumed to be f actor izable into parton distributions in a nucleon ( /tv/a) and the parton 
shadowing factor (Sa/A), 



fa/A{x,Q'^,r) = Sa/A{x,r)fN/A 



(8) 



For the shadowing factor Sa/A we take the parametrization used in the regular HIJING 



model 



30 



55| . We also take into account the intrinsic transverse momentum and the 



transverse momentum broadening due to initial multiple scattering. 

A saturation in the final state, i.e., of p roduced gluons, is possible, even without require- 



ment of saturation in the initial state 



56| . The largest uncertainty in minijet cross sections 



is the strong dependence on the minimum transverse momentum scale cutoff, Pq. It was 
shown that an increased value (with cm. energy ^/s) is required by the experimental data 
indicating that the coherent interaction becomes important. This might be taken also as 
an evidence of gluon final state saturation at very small x [9], lo|. With the Duke-Owens 



parametrization [57|J of parton distribution functions, an energy-independent cutoff scale Pq 
= 2 GeV/c and a constant soft parton cross section a soft = 57 mb are found to reproduce 
the experimental data on the hadron central rapidity density in p + p{p) collisions 36 |. 

However, with a constant momentum cutoff po = 2 GeV/ c, the total number of minijets 
per unit transverse area for independent multiple jet production in central nucleus-nucleus 
collisions, could exceed the limit 

m 



TAA{h)(y]et ^ Po 



(9) 



where T^Aib) is the overlap function of A + A collisions and tt/pq is the intrinsic transverse 
size of a minijet with transverse momentum pq. Therefore, we consider an energy and nuclear 
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size dependent cutoff po{s,A), to ensure the applicability of the two-component model in 
the HIJING model for A + A collisions. The pseudorapidity distribution of charged particle 
in central nucleus-nucleus collisions at RHIC and LHC energies can well be fitted (see Sec. 
IIIA) if we consider a scaling law of the type C A" ^/s^^. 



Pois, A) = 0.416 ^ GeV/c (10) 

New values for po{s, A) increase from po = 1.5 GeV/ c at the c.m.s. energy ^/s = 0.02 TeV 
up to Po = 4.2 GeV/ c at ^/s = 5.5 TeV. These values obtained for central A + A collisions 
are not expected to be valid for peripheral A + A or p + p collisions. These dependences 



are similar to those used in pQCD+saturation model from Ref. 58 



The main difference 



is a factor of two in the constant C, Chij = 0.416 in comparison with Cesk = 0.208 from 
Ref. jssi. Note that Cgsk = 0.208 result in an overestimate of the charged particle density 
at LHC energies by a factor of approximately two jssl . 

The above limit for incoherent minijet production should also depend on impact param- 



eter (b) [59|. In the present calculations within our model such dependence is not included. 
Instead, in the HIJING model an impact parameter dependence of the gluon shadowing is 
considered in the parametrization of the parton shadowing factor Sa/A |55| . Due to shadow- 
ing effects the observed A exponent [a = 0.128) is somewhat less than the number expected 



58 



in the saturated scaling limit {pq{s,A) ~ A^^^) 

The ALICE experiment at the LHC published first experimental data on the charged 
hadron multiplicity density at mid-rapidity in central (0-5%) Pb -|- Pb collisions at 
a/'Snn=2.76 TeV 17|. In this experiment the collaboration also confirmed the presence 
of jet quenching [lot by studying Raa- A measure of the jet quenching at large transverse 
momentum (pt) is given by the ratio of particle yield in A -|- A collisions to that in p + p 
colhsions, defined by 



[l/N^^f)d'NAA/d'pTdy 
N,,n{l/NZ)d^Npp/d^PTdy 

where, A'^evt is the number of events and A'coii is the number of binary nucleon-nucleon 
(NN) collisions. These results provide stringent constraints on the theoretical uncertainties 
in the bulk hadron production and have important consequences on theoretical predictions 
for jet quenching in Pb -|- Pb collisions at LHC energies. 

For a parton a, the energy loss per unit distance can be expressed as dEa/dx = ea/Xa, 
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where ea is the radiative energy loss per scattering and Aa the mean free path (mfp) of the 
inelastic scattering. For_a quark jet at the top RHIC energy [dEq/ dx)Yimc ~ 1 GeV/fm and 



mfp (Ag)RHic ~ 2 fm 37|. The initial parton density is proportional to the final hadron 



multiplicity density. The charged particle density at mid-pseudorapidity at the LHC is a 



factor of 2.2 higher than that at the top RHIC energy [17|. Therefore, for a quark jet at 
the LHC the energy loss (mfp) should increase (decrease) by a factor of ~ 2.0 and become 
[dEql dx)uic ~ 2 GeV/fm and mfp (Ag)LHC ~ 1 fm. Throughout this analysis we will 
consider the results with two sets of parameters for hard interactions: 
{dE/dx)i, i.e., K = 1.5; dEg/dx = 1 GeV/fm; A^ = 2 fm and 
{dE/dx)2, i.e., K = 1.5; dEjdx = 2 GeV/fm; A.^ = 1 fm. 

Since there is always a coronal region with an average length of Xg in the system where 
the produced parton jets will escape without scattering or energy loss, the suppression factor 
can never be infinitely small. For the same reason, the suppression factor also depends on 
Xq. It is difficult to extract information on both dEg/dx and Xq simultaneously from the 
measured spectra in a model independent way 0]. We show below that a constant radiative 
energy loss mechanism (dE/dx=const) as implemented in HIJING/BB v2.0 model provides 
a good description of suppression at intermediate and larger px {4 < px < 20 GeV/ c) for 
charged particle in Pb+Pb collisions at the LHC. 

III. RESULTS AND DISCUSSION 

A. Multiplicity and centrality dependence 

Charged-particle multiplicity density {2dNc\i/ drf) / N^gj-t) in central (0-5%) Pb+Pb colli- 
sions at a/snn=2.76 TeV reported by the ALICE Collaboration is a factor of 2.2 higher 
than that observed in central Au+Au collisions at top RHIC energy. This value is larger 
than most of theoretical predictions, especially those from the CGC models [l^]. Such large 
hadron multiplicity will have iraportant consequences on the estimated shadowing and jet 



quenching in Pb+Pb collisions |63|. The value of {2dNch/ dri) / Np^^t at mid-pseudorapidity 



was reproduced by the HIJINGv2.0 model with a new set of PDFs 9|] and has helped to 



constrain the gluon shadowing parameter in the model 



llQl 



However, these calculations do 



not include the effects of jet quenching, which has been reported recently by ALICE [l9|, 
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ATLAS [21[ and CMS [22 1. 

In previous studies of heavy-ion collisions at RHIC [36|, |37[, the parameters of the HI- 
JING/BB v2.0 model were adjusted in order to reproduce the measured charged particle 
multiplicity as well as (multi) strange particle production. In the present study of heavy-ion 
collisions at the LHC we have to modify the values of these parameters in order to fit ALICE 
data As shown in Fig. [T]a good description (solid curve) of the measured charged par- 
ticle multiplicity density at mid-pseudorapidity is obtained if the parameters po{s, A) = 3.7 
GeV/c, k{s,A) = 6.0 GeV/fm, and {dE/dx)2 (i.e., K = 1.5; dEjdx = 2 GeV/fm; Ag = 1.0 
fm) are used. 
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FIG. 1: (Color online) Comparison of HIJING/BB v2.0 predictions for {2dNch) / dr]) / Np^rt in central 



(0-5%) Pb+Pb collisions at ^snn = 2.76 TeV. The results which give the best fit of the experimental 
value {itI are plotted as a solid curve. The dotted and dashed curves illustrate the sensitivity to 
the main model parameters (see Sec. IIIA for explanation). The dot-dashed curve is obtained for 
p + p NSD collisions at -^/snn = 2.76 TeV. The full triangle is the result obtained by interpolating 



data between 2.36 and 7 TeV in Ref. 



18t |. Comparison with other model calculations are also 



n 



shown: HIJING v2.0 model (open triangle from Ref. [lO|), CGC models (open cross from Ref. 6|]; 
open diamond from Ref. J]). 
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The sensitivity of the calculations to the main parameters is illustrated by the dashed, 
dotted, and solid curves. The dashed curve is obtained with the parameters used at RHIC 
energies: po = 2.0 GeV/c, k,{s) = 2.6 GeV/fm and {dE/dx)i (i.e., K = 1.5; dEq/dx = 1 
GeV/fm; = 2.0 fm). The dotted curve is obtained if po is increased to Pq{s,A) = 3.7 
GeV/ c. Changing from set {dE/dx)i to {dE/dx)2 results in a good description of data (solid 
curve). The dot-dashed curve is the results for p+p non-single diffractive (NSD) collisions at 
y^sjviv=2.76 TeV and are obtained with the same set of parameters as described in Ref. 



The study of the dependence of observable {2dNc\i/ drf) / N^s^^t on the colliding system, 
center-of-mass energy, and collision geometry is important to understand the relative contri- 
butions to particle production of hard scattering and soft processes. Multiparticle produc- 
;ion of charged particles at RHIC energies exhibit the phenomenon of limiting fragmentation 
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40|. This phenomenon is a consequence of Feynman scaling (if one views the collision 



in the rest frame of one of the incident particles, the production process of the soft particles 
is independent of the energy or rapidity of the other particle.). The region in rapidity over 
which the particles appears to be independent of energy, increases with energy. It was ar- 
gued that this phenomenon (called extended longitudinal scaling), is a direct manifestation 



40|. 



of initial saturation which is assumed in CGC models 

Figure 2 display the results for central (0-6%) Au+Au collisions at RHIC energies 
(0.062 TeV < ^/s^ < 0.20 TeV) and in central (0-5%) Pb+Pb coUisions at LHC ener- 
gies (a/snn=2.76 TeV and a/snn=5.5 TeV). Our calculations show that at RHIC energies 
the model predicts approximately a scaling seen also in data. However, some degree of viola- 
tion of the phenomenon of limiting fragmentation and of the extended longitudinal scaling is 
predicted at both LHC energies. This violation is due to the parton hard scattering included 
in our model and missing in CGC models ,|43|. 



Figure 3 shows the dependence of the pseudorapidity densities {{2dNch/ drj) / N^i^^t) at mid- 
pseudorapidity, on the number of participants as measured by the ALICE Collaboration. 
The charged-particle density normalized per participant pair increases by a factor of roughly 
two from peripheral (70-80%) to central (0-5%) Pb+Pb collisions. The model calculations 



which have been tuned to high-energy p + p 



36] and central (0-5%) Pb+Pb collisions data 



reasonably describe the central Pb+Pb data for A'part > 250. However, with the same set of 
parameters the model predicts a stronger rise with centrality than observed. 

For peripheral collisions a better description is obtained if one use as parameters for the 
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FIG. 2: (Color online) Comparison of HIJING/BB v2.0 predictions for central charged particle 
density per participant pair produced in A + A collisions at various energies. The data and the- 
oretical results are plotted in the rest frame of one of the nuclei. The Au+Au data are from the 

The statistical error bars are plotted only at mid-pseudorapidity 



PHOBOS Collaboration 



39 



:br clarity. The Pb+Pb data at -y/sNN=2.76 TeV (filled star) are from the ALICE Collaboration 
l?! . The value at ^snn= 5.5 TeV (star) is obtained by a power-law extrapolation from RHIC 
energies and 2.76 TeV data 

soft part k{s) (instead of k{s,A)), successful in describing p + p collisions, and for hard 
scattering pq{s,A), {dE/dx)i, as deduced from RHIC data. However, in this scenario the 
model overpredicts the absolute magnitude for central collisions (A'part > 250). These results 
indicate that the average energy loss per unit distance and the effective string tension values 
should have an impact parameter (b) and a parton density dependence, and a fine-tuning is 
required to describe better these data. We leave this study to a future analysis. 
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FIG. 3: (Color online) Comparison of HIJING/BB v2.0 model calculations for charged particle 
multiplicity density at mid-pseudorapidity per participant pair {2dNch / dr]) / Npi^j-t as a function of 
number of participants in Pb+Pb collisions at y^SNN = 2.76 TeV. The solid and dotted curves are 
the results obtained using the same set of parameters as in Fig. [1] (see text for explanations). All 
calculations include SCF effects and JJ loops. The data are from ALICE Collaboration 1^ and 
the error bars include statistical and systematic uncertainties. 



B. Nuclear modification factor for charged particle 



In this paper, we make a simultaneous analysis of the high and low pt part of charged- 
hadron pt spectrum. Experiments at RHIC energies reported a suppression by a factor ^ 
5 compared with expectations from an independent superposition of nucleon-nucleon (NN) 



collisions 



61 



62|, interpreted as strong evidence of formation of a strong-coupling QGP 



(sQGP) in A + A collisions. T 
Pb + Pb collisions at the LHC 



le jet 



19 



21 



quenching has also been confirmed in central (0-5%) 



22|. 



Having a good description for charged particle multiplicity at mid-pseudorapidity, we 
analyze the model predictions on suppression of single hadron spectra in central (0-5%) 
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Pb+Pb collisions at the LHC For other recent works discussing NMF of charged particle 
{Raa) at the LHC, see recent references 63N69|. At the higher LHC energies a larger density 
of the medium is expected. This should lead to a larger suppression compared with that 



seen at RHIC energies. The observed nuclear modification factor -Raa |19| is characterized 
by the following behavior: a maximum at approximately px = 2 GeV/ c and a decrease with 
increasing pt in the range 2-6 GeV/c. For pt greater that 6 GeV/c, -Raa rises to reach a 
value of ^ 0.4 at pr ^ 20 GeV/c. 



(a) BBv2,0, 2.76 TeV, h^+h' 

K = 2.6 GeV/fm. (dE/dx)j 

K=6.2GeV/fm, (dE/dx)( 
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(b) BBv2.0. Central 0-5% 

K = 2.6GeV/fm, (dE/dx)^ 

K = 6.2 GeV/fm, (dE/dx)j 

K = 6.2 GeV/fm, (dE/dx)^ 
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FIG. 4: (Color online) (a) Comparison of HIJING/BB v2.0 predictions for the pT distributions 
of primary charged particles at mid-pseudorapidity (—0.8 < r] < 0.8) in central (0-5%) Pb-|-Pb 
collisions at ^/s^ = 2.76 TeV. Theoretical calculations are for different values of parameters (see 
text for explanations). The p + p references values are shown as the lower solid histogram, (b) 
Theoretical calculations for nuclear modification factors -Raa are compared to data. The data (left 
and right panels) are from the ALICE Collaboration {l9]. Error bars are statistical only . In part 
(b) the open triangles are the results for -Raa at pr > 6.5 GeV/c using alternative p + p reference 



spectrum 



The measured and predicted transverse momentum spectra at mid-pseudorapidity of 
primary charged particles in central (0-5%) Pb-|-Pb collisions (upper solid histogram) and 
in p + p collisions (lower solid histogram) at a/snn=2.76 TeV are compared in Fig. Hl^a). 
The HIJING/BB v2.0 model describe fairly well the experimental pt spectrum reported in 
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Ref. 



19| . The sensitivity of the calculations to the effective value of the string tension (k) 



is shown by comparing the dotted = 2.6 GeV/fm) and the dashed ( k{s,A) = 6.2 

GeV/fm) histograms. The sensitivity to the main hard scattering parameters is seen by 
comparing the results obtained with the two sets {dE/dx)2 (solid histograms) and {dE/dx)i 
( dash-dotted histograms). However, with the best set of parameters (solid histograms) 
the model still overestimate the high pt tail of the pt spectra. There are no experimental 
data on charged hadron spectra in p + p collisions at ^/s^ = 2.76 TeV. ALICE data on 
the suppression factor (Raa) were obtained using a. p + p spectrum interpolated from their 
data at 0.9 TeV and 7 TeV. The open triangles in Fig. Sl^b) represent the uncertainty in the 
interpolation method (see Ref. 19| for more details). 

In medium modification of NMFs -Raa is mainly caused by the energy loss suffered by 
partons while transversing the plasma due to collisions and radiation of gluons before they 
fragment and by the nuclear shadowing. To describe the nuclear shadowing we use in Eq. (8) 
the same functions as those used in regular HIJING model 30|. This has been shown to 

. The benchmark 



be successful in explaining the results for Raa at RHIC energies [37], 
of any theoretical calculation is to describe the shape seen at low px < Q GeV/c in -Raa- 
We can explain the shape and the presence of a maximum at approximately = 2 GeV/c 
as a specific interplay between nonperturbative (SCF effects) and perturbative (gluonic JJ 
dynamical loops) mechanism contributions. The results without JJ loops (dashed histogram) 
shown in Fig. |4](b) strongly underestimate the experimental values of -Raa in this pt range. 

Since we are mostly interested in the overall effects, we concentrate also on the modifica- 
tion of high Pt hadron spectra due to an assumed total energy loss related to the averaged 
energy loss per unit distance. The mechanism which gives a constant energy loss per colli- 
sions is seen to describe quite well the data (within the systematics uncertainties) over the 
largest px window of 6 - 20 GeV/c. -Raa increases with pj- because of the constant energy 
loss assumed here. For a constant energy loss per collisions (AEt) the ratio AEt/Et be- 
comes smaller for larger Et, and thus the suppression is expected to decrease with increasing 
Pt- Here, the energy loss per collisions is about twice as large and the mean-freepath of 
parton-parton interactions is smaller by a factor of roughly two as compared with that used 
at the top RHIC energy. These analyses can provide information about the average total 
energy loss the parton suffers during its interaction with the medium. However, a more 
quantitative analysis should be performed with the knowledge of the dynamical evolution 
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of the system which is beyond the scope of this paper. 



C. Baryon/meson anomaly at the LHC 



Identified particles at high pt provide direct sensitivity to differences between quark 
and gluon fragmentation 70|. Proton and pion production at high px is expected to have 
significant contributions from quark fragmentation while antiprotons result mainly from 
gluon fragmentation 60|. Therefore, the ratiosp/vr", p/vr"*" are sensitive to the possible color 
charge dependence of energy loss. Systematic measurements of identified particle spectra in 
p + p and Au+Au collisions at RHIC energies 7CH74l| in the pt region 2-12 GeV/ c show a 



significant suppression with respect to binary scaling. However, protons and anti-protons 
are less suppressed than pions at intermediate pr (2-9 GeV/ c). For pt > 2 GeV/ c, the ratios 
p/7T~ and p/tt~^ do not depend strongly on energy and they are higher in central than in 
peripheral collisions. In addition to being sensitive to quark and gluon jet production these 
are also sensitive to baryon transport properties and the energy density. 

Moreover, the meson/baryon ratio suggest that baryon and antibaryon production may 
dominate the moderate high pr flavor yields. These findings challenge various models incor- 
porating jet quenching 60|, 75|, and/or constituent quark coalescence The fact that 



in the intermediate px region the p/7T~^ and p/n ratios are close to unity ( baryon/meson 
anomaly) has been attributed to either quark coalescence or to novel baryon transport dy- 



namics 



75|, based on topological gluon field configurations |78l-l8l|. However, the baryon 



72|. 



junction model predictions from Ref. [75| are not in agreement with data 

Figure [5] displays the mid-rapidity ratios for central (0-12%) Au-|-Au collisions [72| at 
i/snn = 0.2 TeV in comparison with model predictions. The observed ratios peak at pt ~ 
2 — 3 GeV / c with value close to unity and then decrease with increasing px and approach the 
ratios observed in p + p collisions at > 6 GeV/ c. The results presented (solid histograms) 
include SCF and JJ loops contributions. The calculations based on quark coalescence (dot- 
dashed curves) and thermal models (dotted curves) from Ref. [76| are also plotted. The 



grouping of particle production according to the number of constituent quarks has been 



attributed to quark coalescence from a collective partonic medium 77|. The coalescence 



models 



76 



771] can qualitatively describe these data at intermediate pr < 5 GeV/c, but 



strongly underpredict the results at high pt, due to fragmentation functions used within the 
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FIG. 5: (Color online) Comparison of HIJING/BB v2.0 predictions (solid histogram) for ratios 
p/vr" (left panel) and p/vr"*" (right panel) to STAR data j3] and to coalescence (dot-dashed curves) 
and thermal (dotted cm'ves) model calculations from Ref. [t^. The error bars include only statis- 
tical uncertainties and systematics are introduced by thin continuous curves. 



model. To the contrary, the thermal model calculations [76| strongly overestimate the data 
at high pt- Our model provides an alternative explanation to the baryon/meson anomaly. 
The reason is a specific interplay of contributions from SCF effects and dynamical JJ loops. 
The predictions of our model are in good agreement with data up to the highest pr, in 
the limit of systematics uncertainties. However, the model overestimate the p/n'^ ratio for 
3 < Pt < 6 GeV/c, while it provides a good description for protons (p), it underpredicts 
the pion (vr"'") spectrum in this region. Since the ratio p/n~ is well described by the model, 
these calculations could indicate that the fragmentation functions for production are not 
well estimated. 

With the best set of parameters used to describe charged particle production in central 
(0-5%) and peripheral (70-80%) Pb+Pb collisions at ^/s^ = 2.76 TeV (see Sec. IIIA), 
we present in Fig. [6] the model predictions for centrality dependence of the nonstrange 
baryon/meson ratio (left panel) and of the strange baryon/meson ratio (right panel). Our 
model predicts that these ratios do not depend strongly on energy and centrality. The model 
also predicts that the baryon/meson anomaly will persist up to high p^- For central (0-5%) 
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Pb+Pb collisions at a/snn = 2.76 TeV, in a scenario with jet quenching and shadowing as 
explained in Sec. II, but without J J loops and SCF effects the p/vr" ratio in the region of 
the maximum (2 < < 3 GeV/c) is 0.07. We have studied the sensitivity of the results to 
the JJ loops and SCF effects. Considering only JJ loops leads to an increase by a factor of 
two, and taking into consideration only SCF effects leads to an increase up to a factor of ~ 
14 for the p/vr~ ratio. Note that an important fraction of this increase is also due to much 
stronger quench of pions in comparison with that of antiprotons, as we will show below (see 
Fig. ED. 
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FIG. 6: (Color online) Predictions for centrality dependence of (non)strange baryon/meson ratios 
from HIJING/BB v2.0. The results include SCF effects and J J loops. In each figure the histograms 
correspond to two value of centrality and the associated string tension. 



The ALICE and CMS detectors are designed to perform measurements in the high- 



multiplicity environment of central Pb+Pb collisions at the LHC [3J], 35|. To isolate QGP 
signatures in heavy-ion collisions, understanding the particle production properties from 
A^A^ interactions is necessary. In addition, the transverse momentum spectrum in p + p 
collisions serves as a baseline to study possible suppression (i.e., nuclear modification factor 
-^aaIpt)) in heavy- ion collisions. This stresses the need for reference p + p measurements 
at LHC energies of identified particles. In a previous paper 36| we show that the large 



(strange)baryon-to-meson ratios measured at Tevatron (1.8 and 1.96 TeV) and LHC energies 
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(0.9 and 7 TeV) are well described in the framework of the HIJING/BB v2.0 model. Recently, 
the CDF reported a set of measurements of inclusive invariant pt differential cross sections 
of hyperons with different strangeness, A° (quark content (uds)), S~ (dss), and Q (sss). 
These data are obtained for the central region with pseudorapidity range —l<r]<l and 
Pt up to 10 GeV/ c, and are used to test our model calculations which include SCF effects 
and JJ loops. 

In a scenario with SCF effects (i.e., considering k = k{s) 36|]) the increase of the yield is 



higher for multistrange hyperons (S~,r2~) than for the strange hyperon ( A°), in comparison 
with the results obtained without SCF effects (i.e., taking vacuum string tension value kq = 1 
GeV/fm). The increase of the yield with strangeness content is due to an increase of strange 
quark suppression (7^), obtained from a power-law dependence of effective string tension 
values, n{s) = kq (s/so)°'°^ GeV/fm (see Fig. 2 from Ref. |36|). 
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FIG. 7: (Color online) (a) Comparison of HIJING/BB v2.0 predictions for the Pt differential 
cross sections of A*^ (solid), H (dashed), and $7 (dotted histograms). The calculations are for p + p 



■=r/AO 



collisions at 1.96 TeV and include SCF effects and JJ loops as in Ref. [3g]. (b) The ratios H 
(dashed histogram) and O/A" (dotted histogram). The data are from the CDF Collaboration 
Only statistical uncertainties are shown. 

Figure [7] shows the predicted pt differential cross sections (left panel) for the three hyperon 
resonances in comparison with data. The plots on the right side of Fig. [7] show the ratio 
of the Pt differential cross sections for 'E~ / hP and VL^ / hP. In the limit of the experimental 
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error (systematics+statistical) our model describes well these data. 

Strangeness enhancement, strong baryon transport, and increase of intrinsic transverse 
momentum are all expected consequences of longitudinal SCF effects. These are modeled 
in our microscopic models as an increase of the effective string tension that controls quark- 
antiquark (qq) and diquark-antidiquark (qqqq) pair creation rates and the strangeness sup- 
pression factors (7s). The predictions in central Pb+Pb collisions for initial energy density 
and temperature are eLHC ~ 200 GeV/fm^ and Tlhc ~ 500 MeV, respectively [82|. Both 
values would lead in the framework of our model to an estimated increase of the average 
value of string tension to klhc ~ 5 — 6 GeV/fm {2I, consistent with our parametrization 
from Eq. (5). 
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FIG. 8: (Color online) Predictions from HIJING/BB v2.0 model for nuclear modification factors 
-^AA identified particles at mid-pseudorapidity (—0.8 < ?? < 0.8) in central (0-5%) Pb-|-Pb 
collisions at -^snn = 2.76 TeV. The results include SCF effects and JJ loops. Left: pions, kaons 
and protons. Right: hyperons A, H, fi. 

The predictions for NMFs of identified particles (-RaaIpt)) are presented in Fig. [8] for 
nonstrange particles (left panel) and (multi) strange particles (right panel). From the results 
of our model we conclude that the baryon/meson anomaly, i.e., different meson and baryon 
suppression, will persist in central Pb-|-Pb collisions at the LHC energy ^snn = 2.76 TeV 
up to a higher value of transverse momentum. The -Raa(pt) also exhibit an ordering with 
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strangeness content at low and intermediate pt, the increase of the yield being higher for 
multistrange hyperons than that for (non)strange baryons (i.e., -RpbPb(fi) > -RpbPb(S) > 
-RpbPb(A) > -RpbPb(p))- This prediction could be verified by the future planned experiments 
at the LHC. 



IV. SUMMARY AND CONCLUSIONS 

We have studied the influence of strong longitudinal color fields and of possible multigluon 
J J loops dynamics in particle production in Pb+Pb collisions at 2.76 TeV. We have investi- 
gated a set of observables sensitive to the dynamics of the collisions, covering both longitu- 
dinal and transverse degree of freedom. A detailed comparison with available experimental 
data from the LHC has been performed. 

We find that the inclusion of the multiple minijet source limits the growth of the string 
tension k{s) to be approximately linear as a function of saturation scale Qsat, in contrast to 



recent approaches 43| where k{s) scales as Qsat- The reason is that in the CGC model the 
coUinear factorized minijet mechanism is suppressed by geometric scaling to much higher 
Pt- Moreover, an increase with A due to an increase in initial energy density is necessary in 
order to better describe new ALICE data. 

We have shown that strong color field (SCF) could play an important role in particle 
production at mid-rapidity in heavy-ion collisions. The mechanisms of hadron production 
are very sensitive to the early phase of the collisions, when fluctuations of the color field 
strength are highest. Strong color field effects are modeled by varying the effective string 
tension that controls the quark- ant iquark (qq) and diquark-antidiquark (qqqq) pair creation 
rates and strangeness suppression factors (7^). SCF also modify the fragmentation processes 
resulting in an increase of (strange)baryons. 

We show that the baryon/meson anomaly manifests in an increase of baryon-to- meson 
ratios as well as in different suppression of mesons and baryons and persists up to high p^ 
{ Pt < W GeV/c). We show also that the -RaaIpt) exhibit an ordering with strangeness 
content at low and intermediate pt, the increase of the yield being higher for multistrange 
hyperons than that for (non) strange baryons. 
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